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Abstract. Many problems in computational geometry are not stated in graph- 
theoretic terms, but can be solved efficiently by constructing an auxiliary graph 
and performing a graph-theoretic algorithm on it. Often, the efficiency of the al- 
gorithm depends on the special properties of the graph constructed in this way. We 
survey the art gallery problem, partition into rectangles, minimum-diameter clus- 
tering, rectilinear cartogram construction, mesh stripification, angle optimization 
in tilings, and metric embedding from this perspective. 

1 Introduction 

Graph algorithms and computational geometry form separate communities with sepa- 
rate conferences such as the International Workshop on Graph-Theoretic Concepts in 
Computer Science and the ACM Symposium on Computational Geometry, respectively, 
but they also meet in broader algorithms conferences, and there has been much interplay 
between the research topics in the two areas. 

Many classical graph algorithm problems have geometric analogues, algorithmic 
problems on graphs defined by a geometric input. In most cases, problems of this type 
can be solved directly by constructing the graph and then applying a general-purpose 
graph algorithm, but can be sped up by examining the graph algorithm's structure more 
closely and applying appropriate geometric data structures. A notable instance of this 
phenomenon is the Euclidean minimum spanning tree (the spanning tree of a complete 
graph in which the vertices represent points and edge lengths are Euclidean distances): 
by using a Delaunay triangulation in place of a complete graph, the quadratic time of 
a naive algorithm can be improved to O(nlogn) [67]. Other work along these fines 
includes algorithms for Euclidean matching [1,71] and bipartiteness testing [26]. 

Graph drawing, the visualization of graphs via geometric graph representations [19, 
45, 59], forms another community represented by the annual International Symposium 
on Graph Drawing. Most work in the area concerns drawings in which a graph's ver- 
tices are represented as geometric points, disks, or polygons, and its edges are repre- 
sented as fine segments or curves. Researchers in this area seek algorithms that optimize 
mathematical stand-ins for their aesthetic quality and readability such as the number of 
crossings, the number of bends in non-straight edges, the angles formed by adjacent 
edges, the area of the drawing, or the amount of syiimietry that the drawing displays; 
the interplay between these different measures provides much scope for research. 

Geometric techniques have sometimes also been applied to solve problems that were 
originaUy defined in purely graph-theoretic terms. One instance concerns parametric 



minimum spanning trees, minimum spanning trees of graphs whose edge weights are 
linear functions of a parameter. Dey proved [18] that 0{mn^^^) different trees may be 
formed in this way from a parametric graph with m edges and n vertices, using the 
crossing number inequality that s simple plane curves with n shared endpoints have 
Q.{s^ /n^) crossings [4,50], while the best known lower bound, Q.{ma{n)), involves a 
reduction from another geometric problem, lower envelopes of line segments [22]. 

In this paper we survey connections between computational geometry and graph 
algorithms of yet another type: problems in computational geometry that, although not 
initially defined graph-theoretically, may be solved by constructing an auxiliary graph 
from the input, applying a purely graph-theoretic algorithm to this auxiliary graph, and 
translating the output of this algorithm back into geometric terms. We discuss prob- 
lems including the art gallery problem, partition into rectangles, minimum-diameter 
clustering, rectilinear cartogram construction, mesh stripification, angle optimization in 
tilings, and metric embedding. The ordering of the problems is roughly chronological, 
and the selection of topics is idiosyncratic and (especially for the later problems) largely 
drawn from the author's own research rather than being comprehensive. 

2 Art gallery problems 

Most computer scientists are familiar both with Chvatal's art gallery theorem [13, 61] 
that every «-vertex simple polygon has a set of [n/3j guard points from which the 
whole polygon may be seen, and with Fisk's elegant graph-theoretic proof [36]. One 
begins the proof by adding a maximal set of non-crossing diagonals to the polygon, 
partitioning it into triangles. Graph-theoretically, the vertices, sides, and added diago- 
nals of the polygon form a maximal outerplanar graph; the weak dual of this graph (the 
adjacency graph of the triangles, omitting the outer face) is a tree (Figure [T] center). 
Every maximal outerplanar graph may be colored with three colors, as may be proven 
by induction: the result is clearly true when the graph is a triangle, and any larger max- 
imal planar graph may be colored by removing a leaf from the dual tree, coloring the 
remaining graph, restoring the leaf, and observing that the restored vertex has only two 
neighbors and therefore has a free color available. Each triangle has one vertex of each 
color, so each color class forms a valid guard set, and the smallest of the color classes 
has at most [n/3j vertices. This proof technique translates into an efficient algorithm: 
polygons may be triangulated in linear time [10], and the induction proof leads to a 
linear-time algorithm for 3-coloring maximal planar graphs. 

A less well-known variant of the art gallery problem concerns simple orthogonal 
polygons, simple polygons all of whose sides are parallel to the coordinate axes. Sim- 
ple orthogonal polygons need only [n /4J guards, as can be shown again by graph col- 
oring. Every simple orthogonal polygon can be partitioned by diagonals into convex 
quadrilaterals [46]; the resulting tree of quadrilaterals can be viewed as a special type 
of squaregraph, a planar graph in which every bounded face has four sides and every 
vertex either belongs to the unbounded face or has at least four incident edges [6, 12]. 
Adding the two diagonals of each quadrilateral to a squaregraph forms a kinggraph [12] 
(Figure [T| right). As with maximal outerplanar graphs, the kinggraphs derived from 
simple orthogonal polygons may be shown to be 4-chromatic by removing leaves of the 



Fig. 1. The art gallery problem. Left: a comb polygon requiring [«/3j guards. Center: 
triangulating an input polygon produces a maximal outerplanar graph, the weak planar 
dual of which is a tree. Right: The kinggraph formed by adding diagonals (dashed) to a 
partition of an orthogonal polygon into convex quadrilaterals. 

dual tree|^In any four-coloring, each quadrilateral has a vertex of each color, so each 
color class forms a guard set, and the smallest color class has at most [n/4j guards. 
This bound is tight: a comb-shaped orthogonal polygon requires [n/4j guards, just as 
a comb-shaped simple polygon requiring [n/3j guards shows that the bound for the 
standard art gallery problem is tight (Figure [T] left). 

The biggest remaining open problem in art gallery theory concerns edge guards: 
how many edges must one select, from a simple polygon, so that every point of the 
polygon may be seen from some point on some selected edge? It is not clear, in this 
case, what graph one should define on the edges of the polygon to force the color classes 
of a coloring of the gi'aph to form edge guard sets. 

3 Partition into rectangles 

Many geometric algorithms take as input a complicated polygonal domain and cover 
or partition it using simpler shapes [65]; the partitions into triangles and quadrilaterals 
from the previous section are important examples. Another problem of this type con- 
cerns the partition into rectangles of an orthogonal polygon. The input polygon may 
have polygonal holes, unlike the simple orthogonal polygons of the previous section, 
and the rectangles are not required to meet edge-to-edge and vertex-to-vertex. The goal 
is to minimize the total number of rectangles in the partition (Figure|2]l. 

Rectangular partitions have many applications. In VLSI design, it is necessary to de- 
compose masks into the simpler shapes available in lithographic pattern generators [62], 
and similar mask decomposition problems also arise in DNA microarray design [39]. 
Rectangular partitions can simplify convolution operations in image processing [35] 
and can be used to compress bitmap images [11]. Closely related matrix decomposi- 
tion problems have been applied to radiation therapy planning [21,47], and rectangular 
partitions have also been used to design robot self-assembly sequences [51]. 

Define a good diagonal to be an axis-parallel line segment interior to the input 
polygon that connects two concave vertices of the polygon. As several authors inde- 



More generally, every kinggraph is 4-chromatic, but a proof is beyond the scope of this survey. 




Fig. 2. Partitioning an orthogonal polygon (left) into the minimum number of rectan- 
gles (center). The right figure shows the axis-parallel diagonals that connect pairs of 
concave vertices; the rectangle partition problem may be solved by finding a maximum 
independent set in the bipartite intersection graph of these diagonals. 



pendently discovered [35, 55, 60], the minimum number of rectangles in a partition of 
a polygon with n vertices and h holes is«/2 + /! — g — 1, where g is the maximum size 
of a set of disjoint good diagonals. To see this, consider the number of corners of rect- 
angles in a partition, four times the number of rectangles. Let G be a maximal set of 
disjoint good diagonals that form a subset of the line segments in some partition of an 
orthogonal polygon, and define a bad vertex to be a nonconvex polygon vertex that is 
not an endpoint of G. Each polygon vertex forms at least one rectangle comer. Addi- 
tionally, for each bad vertex v, let s be an interior line segment of the partition having 
V as an endpoint. Then either s crosses a segment in G, and the two corners formed by 
the crossing that are on the same side as v with respect to the crossed segment can be 
charged to v, or s ends at a non- vertex and the two corners formed at its other endpoint 
can be charged to v. This charging scheme shows that the number of rectangle corners 
in the partition is at least n plus twice the number of bad vertices. The number of non- 
convex vertices in any orthogonal polygon with h holes is «/2 + 2/; — 2, so the number 
of bad vertices isn/2 + 2/i — 2|G|— 2, and a lower bound ofn/2 + h—\G\ — \ rectangles 
follows. Conversely, if G is any set of disjoint good diagonals, a partition with exactly 
n/2 + h— \G\ — 1 rectangles may be found by considering the bad vertices for G, one 
at a time, and extending a fine segment from each bad vertex to the closest previously- 
added segment or polygon side. Thus, finding a partition into a minimum number of 
rectangles is equivalent to finding a maximum number of disjoint good diagonals. 

The intersection graph of the good diagonals is bipartite: each horizontal diago- 
nal intersects only vertical diagonals and vice versa. Therefore, finding the maximum 
number of disjoint good diagonals translates, in graph-theoretic terms, into finding a 
maximum independent set in a bipartite graph. By Konig's theorem [49], in any n- 
vertex bipartite graph the maximum independent set has size n — M, where M is the 
cardinality of a maximum matching; an independent set of this size may be found from 
a maximum matching by partitioning the vertices according to the lengths of the short- 
est alternating paths from an unmatched vertex to the given vertex, and including the 
vertices at even levels of this partition. Therefore, the maximum independent set of the 
intersection graph, the corresponding maximimi set of disjoint good diagonals, and a 
partition into a minimum number of rectangles, may all be found in polynomial time. 



A naive implementation of an algorithm that reduces the problem to a bipartite 
graph and then applies a general-purpose bipartite graph matching algorithm would 
solve the problem in time 0{n^'^), where n denotes the number of vertices of the input 
polygon: this is the time needed to apply the Hopcroft-Karp matching algorithm [41] to 
the intersection graph of the good diagonals, which may have 0(n) vertices and 0(n^) 
edges. However, by using geometric range searching data structures to speed up the 
search for alternating paths within the matching algorithm, it is possible to improve the 
overall running time to 0{rr'l'^\ogn) [42,53,54]. 

4 IVIinimum diameter clustering 

The diameter of a set of points is the maximum distance of any two of its points. The 
problem of finding low-diameter subsets of larger sets of input points [3, 27] may be 
formulated in several different ways: one may take as input a number k and produce as 
output the subset of k points with minimum diameter, one may take as input a number 
D and produce as output the largest subset with diameter at most D, or one may solve 
the decision problem in which D and k are both given and the task is to determine 
whether there exists a set of k points with diameter at most D. Since there are only 
O(n^) potential diameter values and n different values of k, an efficient algorithm for 
any one of these tasks leads to efficient algorithms for the other two tasks as well. 

Finding the largest size A: of a set with diameter at most a given value D has a 
direct translation to a graph theoretic problem, but the very directness of the translation 
means that it is unhelpful: it provides merely a restatement of the problem rather than 
conveying any new insight. We may scale the input point set so that D — 2; then what 
we seek is the largest clique in a unit disk graph, the intersection graph of unit disks 
centered at the points (Figure|3] left); note that a set of disks forming a clique need not 
have a common intersection point. Maximum cliques in unit disk graphs may be found 
in polynomial time, given a disk representation of the graph [15] but this is just a trivial 
restatement of the minimum diameter clustering problem; it is NP-hard to find a disk 
representation given only a graph-theoretic description of a unit disk graph [8]. 

However, just as in the rectangle partition problem, minimum diameter clustering 
may be reduced to a more basic graph-theoretic problem, maximum independent sets 
in bipartite graphs [3, 15]. Suppose that we know or can guess the two points p and q 
forming the diametral pair in a cluster. Then, all the other points of the cluster must be 
within the lime formed by intersecting the two circles with pq as radius, centered on p 
and q. The largest possible cluster within the lune is the maximum clique of the unit 
disk graph with disk centers inside this lune, or equivalently the maximum independent 
set of the complement G of this restricted unit disk graph. But G is bipartite: if the lune 
is bisected by line pq, then a point in G can be connected by an edge only to other 
points on the other side of the bisection line (Figure [3] right). Therefore, its maximum 
independent set may be found in polynomial time, as discussed in the previous section. 

Based on this idea, one could test all 0(«^) pairs pq, find the bipartite graph derived 
from each pair and apply a bipartite matching algorithm to it, and return the best cluster 
found from all of these separate tests. However it is more efficient to use dynamic graph 
algorithms to share work between multiple different matching problems [27]. Suppose 



Fig. 3. Minimum diameter clustering. Left: a unit disk graph and its largest clique (the 
five darker circles in the upper left). Right: for any two specified points, the largest 
cluster having those two points as its diametral pair may be found as a maximum inde- 
pendent set of a bipartite graph in which the two sides of the bipartition are two halves 
of a lune defined by the two specified points. 



we seek the size of the largest cluster for a given D. For each point p that could be 
an endpoint of a diametral pair, consider the lune defined by a segment of length D 
with p as one endpoint. If the defining segment rotates through an angle of 2% aroimd 
p, the lune rotates with it, and the set of input points inside the lune undergoes a se- 
quence of 0{n) discrete changes in which some point joins or leaves the set. After each 
change, one may update the maximum matching of the bipartite graph defined from 
the lune by a single alternating path search. Thus, the overall algorithm loops through 
all n possible choices of p, and performs a nested loop through the 0{n) set insertion 
and deletion operations defined by rotating a size-D lune around p. For each set up- 
date operation the algorithm performs an alternating path search to update a maximum 
matching and the maximum independent set in the bipartite graph defined from the lune, 
and when the nested loops terminate the algorithm returns the largest cluster found in 
each of these searches. As Aggarwal et al. [3] describe, each step of an alternating path 
search may be performed in logarithmic time with the aid of the circular hull data struc- 
ture of Hershberger and Suri [40]. Therefore, each alternating path search takes time 
O(nlogn), the sequence of 0{n) altemating path searches for a single choice of p takes 
time O(n^logn), and the overall clustering algorithm takes time 0{ir'\ogn) [27]. 

If k rather than D is given as input, the problem may be solved by a binary search 
among the (9(n^) different distances defined by the input points, that checks for each 
distance whether it is the diameter of some A:-point cluster. For this variant, the time may 
be further improved to 0{n logn -|- k^nlog^ k) by using a fe-nearest-neighbor calculation 
to reduce the problem to 0{n/k) subproblems with 0{k) points per subproblem [27]. 

It appears to be open whether there exist polynomial time algorithms for solving the 
analogous minimum diameter clustering problems in higher dimensions (equivalently, 
finding maximum cliques in intersection graphs of unit balls in higher dimensions) or 
for finding a maximum clique in a unit disk graph when a geometric representation of 
the graph is not given. 



Fig. 4. A rectilinear cartogram of five countries in northwestern Eu- 
rope. Map based on a CC-BY-SA licensed image by Brianski, 
Canuckguy, Zaparojdik, Madman2001, Roke, and Ssolbergj, online at 
http://commons.wikimedia.Org/wiki/File:BlankJMap-ofJEurope_-w_boundaries.svg 



5 Bend minimization 

A rectilinear cartogram [58,64,72] is a diagram in which geographic regions have been 
replaced by orthogonal polygons, with approximately the same shapes and in approx- 
imately the same positions with respect to each other as they hold geographically, but 
in which the areas of the regions have been distorted to represent numerical data about 
the regions unrelated to their physical areas. In introducing these diagrams in 1934, 
Raisz [64] wrote, "it should be emphasized that the statistical cartogram is not a map," 
and the stylization inherent in using orthogonal polygons helps perform this emphasis. 

We have studied algorithms for constructing cartograms that can accomodate arbi- 
trary area assignments [31] and in which the adjacencies between regions have desired 
orientations [30] but a more basic problem, constructing a cartogram in which the re- 
gion shapes are as simple as possible, had already been solved in a different context, that 
of graph drawing [69]. Simple shapes such as rectangles aid the viewer in measuring 
and comparing areas and hence in understanding the data represented by a cartogram. 
We would like to find a cartogram that represents a given map with a minimum number 
of bends (corners interior to the boundary between two adjacent regions); for instance, 
in Figure |4] there is one bend, between Belgium and Luxembourg; this is optimal, as 
Luxembourg is surrounded by only three countries and hence must have at least one 
bend on its boundary. As Tamassia [69] shows, this problem of bend minimization can 
be solved in polynomial time by translating it into a network flow problem. 

To form a flow problem that represents the bends of some given cartogram, construct 
a graph that has a single "circulation" vertex, a vertex for each region of the cartogram 
(including the exterior of the diagram as one of the regions), and a vertex for each 
junction where three or four regions of the cartogram meet (four regions at a junction 
may be needed to model places like the Four Comers in the southwest of the U.S. 
where four states meet, but five or more regions cannot meet at a single junction in an 
orthogonal cartogram). Include edges from the circulation vertex to each other vertex, 
between each two vertices that represent adjacent regions, and between two vertices 
that represent an incident junction-region pair. Label this graph with flow amounts. 



as follows: each junction vertex has four incoming units of flow from the circulation 
vertex, and sends one unit of flow to the vertices representing the regions in the four 
quadrants incident to the junction. For each bend in the cartogram, send one unit of flow 
from the region that has a convex vertex at the bend to the region that has a concave 
vertex. For each interior region having k junctions on its boundary, send A — 2k units 
of flow to the circulation vertex (or equivalently, if 4 — 2A: < 0, send 2k — A units from 
the circulation vertex to the region vertex). And for an exterior region with k junctions 
on its boundary, send 2k + 4 units of flow to the circulation vertex. The result can be 
shown to be a vaUd circulation: that is, for each vertex of the flow graph, the numbers 
of incoming and outgoing flow units are equal. For the junction vertices, this is clear 
because the four incoming units are balanced by the four quadrants into which a unit 
of flow is sent. An interior region with A: junctions and no bends must form a rectangle, 
as a junction cannot form a concave comer; thus, it has four incident junction vertices 
that send a single unit of flow, and the remaining k — 4 incident junctions send two units 
of flow into the region, for a total of 2A: — 4 incoming units, balancing the flow to the 
circulation vertex. Each additional cave comer at a bend causes one unit of incoming 
flow from another region, but must be balanced by an additional convex corner; if that 
convex comer belongs to a bend, it provides a unit of outgoing flow, and if it belongs 
to a junction then that junction sends one fewer unit of incoming flow, in either case 
leading to the same total flow balance. A similar argument shows that the flow into and 
out of the exterior region vertex is also balanced, from which it follows that the flow 
must be balanced at the single remaining vertex, the circulation vertex. 

Conversely, as Tamassia shows, one can assign costs and flow capacities to this 
network in order to force a minimum cost circulation to correspond to a minimum-bend 
drawing. Capacity constraints are needed to force the incoming flow to each junction 
vertex to be exactly four units; the edges between adjacent region vertices are assigned 
unit cost, and all other edges have zero cost. With these constraints and costs, the flow 
described above has cost equal to its number of bends. Conversely, any integer solution 
to the capacitated circulation problem can be translated into a drawing in which the 
total number of bends is equal to the cost of the circulation, so a cartogram with the 
minimum number of bends can be found from a minimum-cost circulation, which can 
be constracted in polynomial time. The flow graph is an apex graph: if one vertex, the 
circulation vertex, is removed, the rest is planar. Therefore, specialized techniques for 
finding flows in planar graphs may be used to speed up this algorithm [38, 43, 57]. 

6 Mesh stripification 

Stripification refers to the problem in computer graphics of partitioning a triangulated 
surface model of a three-dimensional object into triangle strips, sequences of triangles 
meeting edge-to-edge [5, 29, 34, 73]. Such a partition allows for fast rendering by re- 
quiring the coordinates of only one vertex per triangle to be transmitted to the graphics 
hardware; the other two vertex locations may be found from the previous triangle in the 
strip. Triangle strips aid in data compression of geometric models: a predicted location 
that aids in compressing the coordinates of each successive vertex may be found by 
extrapolation from the previous triangle in a strip [17]. AdditionaUy, if a mesh can be 



Fig. 5. Left: a single strip in a triangulated torus, rendered by M. Gopi, from [25]. Right: 
a local move at a vertex allows three cycles to be merged into a single cycle. 



represented as a single cyclic triangle strip, as in Figure |5] (left), the structure of the 
strip may be used to cover the surface of the model with a space-filling curve [29], 
and contractions of the boundary edges of the strip can be used for topology-preserving 
simphfications of the model [20] . 

Finding such a strip may be viewed graph-theoreticaUy, as finding a Hamiltonian 
cycle in the 3-regular dual graph of the triangulation, a graph with one vertex per tri- 
angle and one edge for each pair of adjacent triangles. However, this is problematic 
for two reasons. First, the Hamiltonian cycle problem is NP-complete and the known 
exponential-time algorithms for the problem [25, 44] are only capable of solving it 
within a reasonable time for models of at most a few hundred triangles. Second, and 
more importantly, not all triangulated models, even with spherical topology, have sin- 
gle triangle strips of this type. Tutte's counterexample to Tait's conjecture that 3-regular 
polyhedra are Hamiltonian [68,70] dualizes to become a triangulated mesh that cannot 
be represented as a single cycle of triangles. 

However, it is important to realize that the dual graph does not uniquely represent 
the geometry: we may change the triangulation, and hence change the graph, without 
changing the model's shape. The dual graph of any triangulated model is both 3-regular 
and bridgeless; therefore, by a theorem of Petersen [63], it has a perfect matching, which 
may be found efficiently [7]. The set of edges complementary to the matching forms 
a partition of the triangulation into a collection of disjoint cycles [29]; however, there 
may be more than one of these cycles. As we observed [29], in many cases a local move 
at a vertex of the triangulation, that swaps matched and unmatched edges connecting 
the triangles sharing that vertex, may reduce the number of cycles (Figure |5] right). 
If no such move is available, two adjacent triangles from two different cycles may be 
bisected across their shared edge, leading to a new triangulation of the same surface 
with the property that a cycle-reducing local move is available at the newly created 
vertex. By repeating this subdivision process, one eventually reaches a triangulation 
that covers a surface identical to the input model, but one that has a different dual graph 
than the input and that can be represented as a single strip. Although theoretically this 



Fig. 6. Two drawings of the squaregraph dual to a 5 -chromatic triangle-free circle graph 
described by Ageev [2], from [6]. 

could increase the total number of triangles by as much as a factor of 3 /2, in practice 
we saw at most a 2% increase. 

7 Angle optimization of tilings 

The planar dual of any arrangement of lines in the plane may be represented as a tiling 
of a convex subset of the plane by centrally symmetric convex polygons; this tiling 
has one polygon for each crossing point of the arrangement, which may be chosen to 
have unit-length edges perpendicular to the lines crossing at that point. De Bruijn [16] 
used this construction to form aperiodic tilings (including the rhombic Penrose tiling) 
from overlaid families of evenly-spaced parallel lines. More generally the dual of a 
weak pseudoline arrangement (a set of curves that are the image of a line under a 
homeomorphism of the plane, that have at most one point of intersection per pair of 
curves, and that cross each other at their intersection points), may be represented as 
a tiling of a nonconvex simple polygon by centrally-symmetric polygons, and every 
such tiling arises in this way [24,28]; this provides a convenient method for visualizing 
squaregraphs, which are the duals of triangle-free hyperbolic line arrangements [6]. 

However, the tilings obtained directly from this construction may be hard to read 
due to having polygons with very sharp angles, as is true for instance for the tiling 
in Figure |6] (left). In unpublished work with Kevin Wortman [32], we considered the 
problem of finding a combinatorially equivalent tiling that maximizes the minimum 
angle in the tiling (the so-called angular resolution [9,56]), such as the one in Figure|6] 
(right); we showed that this optimal tiling could be found by a parametric shortest path 
computation in an auxiliary graph derived from the input. 

Our algorithm constructs a graph in which the edges have weights that are linear 
functions of a parameter X that represents the angular resolution of the drawing. In a 
tiling of this type, one can define an equivalence relation on the sides of tiles in which 



opposite pairs of sides on the same side are equivalent; the equivalence classes form 
zones of line segments that are required to be parallel and to have the same length. Our 
graph has one vertex per zone, and an additional start vertex that has a zero-length edge 
to each other vertex; the distance from the start vertex to a zone's vertex represents the 
adjustment in angle for the segments in the zone between an initial tiling and the optimal 
tiling. The constraint that the angle between two sides of a polygon be at least A, can then 
be expressed by the existence of an edge between the vertices representing the zones 
containing the two sides, with length 9, — Qj — X, where 9, and Qj are the angles formed 
by the two zones in the initial tiling. Similarly, the constraint that an interior angle be 
convex can be expressed by an edge with length 7t+ 9, — Qj. With these vertices and edge 
lengths, it can be shown that there exists a tiling with angular resolution at least X if and 
only if the result of substituting that value into each edge weight function is a graph with 
no negative cycles. Therefore, the optimal angular resolution is the largest value of X 
giving no negative cycles. Due to the special form of the weights in the parametric graph 
(each weight is either constant or a constant minus X) this parametric negative cycle 
detection problem can be solved in (9(n^) time by an algorithm of Karp and Orlin [48]. 
The translation from the tiling angular resolution problem to the parametric negative 
cycle detection problem and back can be performed within the same time bound. 



8 Metric embedding into stars 

There has been a large amount of interest recently within the theoretical computer sci- 
ence community in problems of embedding unstructured metric spaces (which may be 
specified, for instance, as a distance matrix) into simpler and more highly constrained 
metrics, with low distortion [52]. Such embeddings may be used, for instance, in ap- 
proximation algorithms: one can design an approximation algorithm for the constrained 
class of metrics, and apply it to arbitrary metrics, incurring the distortion of the embed- 
ding as a penalty factor in the approximation ratio of the algorithm. The construction of 
graph spanners [23] may also be seen as an instance of this type of problem: one wishes 
to approximate a metric describing the shortest paths in an arbitrary graph, by a more 
highly constrained metric of shortest paths in a sparse graph. Most work in this area 
has concentrated on finding embeddings that guarantee low but non-optimal distortion; 
however there has also been some work on finding the best possible embedding [33,37]. 
In particular, we describe here our work on finding optimal embeddings into star met- 
rics, which (as in the angle optimization of the previous section) involves translating 
the problem into one of parametric negative cycle detection in an auxiliary graph. 

Following [33], we define a star metric to be a metric space in which there exists a 
distinguished point h (the hub) of the metric) such that, for every two points p and q, h 
lies on a shortest path from p to q. Expressed algebraically: d{p,q) = d(p,h)+d{h,q). 
Thus, in contrast to arbitrary «-point finite metric spaces which have n(n — l)/2 degrees 
of freedom, a star metric has only n—l degrees of freedom. We wish to find a map from 
an arbitrary finite metric space (represented as a distance matrix D) to a star metric space 
(represented as a vector H of distances from each point to the hub) that minimizes the 
distortion of the embedding; we do not require that the hub be the image of any point 
in the input space. By scaling H, if necessary, we can assume without loss of generahty 



Fig. 7. The parametric network into which an optimal star embedding problem is trans- 
lated. From [33]. 

that the map does not decrease any distance; thus, what we seek is a vector H satisfying 
the constraint that, for all p and q, D[p^q] < H[p] +H[q], and minimizing the dilation 

5 = max^^[4±^. 
p.q D\p,q\ 

This is a linear program: it may be rephrased as the problem of finding H and 5 sub- 
ject to the linear nondecreasing-distance constraints and the linear constraints that, for 
each p and q, D[p,q]?) — H[p] —H[q] > 0, and, among all tuples of H and 5 that satisfy 
the constraints, finding the tuple minimizing 5. However, it may be solved more effi- 
ciently, in strongly polynomial time, by transforming it into a parametric negative cycle 
detection problem. 

Specifically, as our paper shows, we may find the optimal star embedding using an 
auxiliary graph defined from the input, having two vertices p and p for each metric 
space point p, as well as a special start vertex s. There is an edge with length zero from 
s to each vertex p and from each vertex p to the corresponding p. In addition, for each 
pair of points p and q, there is an edge of length —D[p,q] from p to q and another edge 
of length XD[p,q] from p to q. As we show, the optimal dilation 5 is the smallest value 
that can be assigned to A, such that the resulting graph has no negative cycles, and for 
this value of X the distance D[p] of p from the hub in the optimal embedding may be 
computed as half of the difference between the two distances from s to p and p. 

Thus, the star metric embedding problem may be solved optimally by a single para- 
metric negative cycle detection calculation similar to that in the previous section. How- 
ever, in this case the edge weight functions of the parametrized network no longer have 
the special form needed to apply the Karp-Orlin algorithm; instead, we developed a 
more general algorithm that solves any parametric negative cycle detection problem in 
strongly polynomial time. The basic idea of the algorithm is to use a matrix squaring 
method of Savage [66] to calculate the parametrized function representing the weights 
of paths between each pair of vertices; after / iterations of squaring a matrix whose en- 



tries represent these functions, they will correctly describe the lengths of the shortest 
path between each pair of vertices among paths constrained to have at most 2' hops. 
After each matrix- squaring step, the entries of the matrix will not be Unear functions 
but rather more general piecewise Unear concave functions; we perform a binary search 
among the set of breakpoints of these functions to narrow down the range of values of 
X within which the optimal value must he, allowing us to simphfy the functions to be 
linear once again. Each binary search step involves applying the Bellman-Ford algo- 
rithm to determine whether a specific value of X causes the parametrized graph to have 
a negative cycle. The negative cycle detection algorithm resulting from this parametric 
search procedure has total nmning time 0{n^ log^ «), which is therefore also the time 
for using this approach to solve the optimal star metric embedding problem. 

9 Conclusions 

As we have seen, a graph-theoretic point of view can be useful in many algorithmic 
problems that do not, initially, seem to have much to do with graphs. The graph prob- 
lems occurring in the solutions of the geometry problems discussed here are mostly 
well-known and classical: they include the maximum independent set and maximum 
clique problems, maximum cardinality matching, shortest paths, and minimum cost 
flow. However, sometimes these problems occur in a somewhat different form than they 
have been studied elsewhere in the literature: the star metric embedding problem, for 
instance, required parametric negative cycle detection at a level of generality that had 
previously not been considered. 

Special classes of graphs, and specialized algorithms that take advantage of these 
special classes, seem to be ubiquitous in this area. In the examples here, we have seen 
natural applications of trees, maximal outerplanar graphs, squaregraphs and kinggraphs, 
bipartite graphs, unit disk graphs, apex graphs, and 3-regular bridgeless graphs. The 
simple coloring algorithm for maximal outerplanar graphs described in the section on 
the art gallery theorem generalizes to chordal graphs (maximal outerplanar graphs are 
chordal) and to perfectly orderable graphs [14]. 

There are undoubtedly many more applications of graph-theoretic concepts in com- 
putational geometry waiting to be discovered, and it is important that the fines of 
communication remain open between the graph algorithm and computational geometry 
communities, so that computational geometers will know where to find the algorithms 
they need and so that graph algorithms researchers can focus their efforts on the prob- 
lems and graph classes with the greatest benefit in geometric appfications. 
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